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摘要 :【 目 的】 筛选 特定 条 件 下 桃 星 Conogethes punctiferails 稳定 表达 的 内 参 基因 , ARERR A 
定量 研究 黄 定 基础 。【 方 法 】 参 考 文献 报道 ,并 根据 桃 星 蝇 转 录 组 测序 结果 筛选 出 B- 肌 动 蛋白 基因 
4C 人 甘油 醛 -3- 磷 酸 脱 氢 酶 基因 GAPDH、 核 糖 体 蛋 白 49 基因 RP49 、a- 微 管 蛋 白 基 因 a-tub、 核 糖 体 
蛋白 L13 基因 RPL13 和 液 泡 型 ATP 合 酶 基因 V-ATPase 等 6 个 候选 基因 。 再 利用 实时 荧光 定量 
PCR(qRT-PCR ) 技 术 测 定 6 个 候选 基因 在 桃 蛙 嵌 不 同 发 育 时 期 和 成 虫 不 同 组 织 中 的 表达 量 ; 然后 
通过 ACt 法 ,GeNorm，NormFinder BestKeeper 和 RefFinder 5 个 软件 对 6 个 候选 基因 的 稳定 性 进行 
评估 ;最 后 以 桃 星 蝇 嗅觉 受 体 基因 (olfactory receptor co-receptor, Orco) 和 性 信息 素 结 合 蛋 和 白 基 因 
(pheromone binding proteinl, ，PBP1 ) 为 目的 基因 ,验证 6 个 候选 基因 分 别 在 桃 星 螨 不 同 发 育 时 期 和 
成 忠 不 同 组 织 中 的 的 稳定 性 。[【 结果 】ACt 法 ，GeNorm NormFinder 和 BestKeeper 分 析 结 果 表 明 ,4 
个 软件 对 6 个 候选 基因 稳定 性 的 排序 相似 ,在 桃 昱 崇 不 同 发 育 时 期 和 成 虫 不同 组 织 中 稳定 性 较 高 
的 3 个 基因 为 RP49，RPLI3 和 CAPDH; 同 时 4 个 软件 均 判 定 ACT 的 稳定 性 最 低 。 进 一 步 利 用 
RefFinder 进行 综合 分 析 , 结果 显示 ,在 桃 星 剖 成 求 不 同 组 织 中 ,最 稳定 的 基因 为 RPL13 ,其 次 为 
RP49; 而 在 不 同 发 育 时 期 ,最 稳定 的 基因 为 RP49 ,其 次 为 CAPDH。 另 外 经 CeNorm 软件 计算 得 出 ， 
在 桃 昱 剖 不 同 发 育 时 期 和 成 虫 不 同 组 织 中 桃 星 蜡 内 参 基因 的 最 佳 数目 为 2。 最 后 ,以 噢 觉 受 体 基 
Orco 和 PBP1 为 目的 基因 ,对 筛选 出 的 候选 基因 的 稳定 性 进行 验证 ,发 现 当 使 用 RPL13 和 RP49 
以 及 RP49 和 GAPDH 作为 内 参 基因 时 ,Orco 及 PBPI1 的 表达 量规 律 一 致 ,并 与 桃 昱 剖 生 活 习 性 及 前 
人 研究 结果 一 致 ;而 使 用 ACT 作为 内 和 参 基因 计算 得 到 的 Orco 和 PBP1 表达 量 则 不 规律 。【 结论 ] 在 
不 同 发 育 时 期 桃 昱 蝇 基 因 定 量 研究 中 建议 以 RP49 和 GAPDH 作为 内 参 基 因 , 而 在 桃 昱 蝇 成 虫 不 同 
组 织 中 基因 定量 研究 中 建议 以 RPL13 和 RP49 作为 内 参 基因 。 

关键 词 : MI, 内 参 基 因 ; EO put. 表达 稳定 性 ; 实时 荧光 定量 PCR 
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Screening of reference genes for qRT-PCR in Conogethes punctiferails 
( Lepidoptera: Crambidae) 

YANG Ling, HU Xiao-Jing, XU Zhi-Feng, HE Lin, XIAO Wei ' (College of Plant Protection, Southwest 
University, Chongqing 400716, China) 

Abstract: [ Aim] This study aims to screen and verify stably expressed genes under given conditions as 
reference genes for quantitative real-time PCR (qRT-PCR) in Conogethes punctiferalis, so as to provide 
the basis for quantitative studies of genes of this moth. [Methods] Based on transcriptomics sequencing 
results in C. punctiferalis and reported reference genes in other insect species, six candidate genes 
including B-actin gene (ACT), glyceraldehyde 3-phosphate dehydrogenase gene (GAPDH) , ribosomal 
protein 49 gene ( RP49) , alpha tubulin gene (a-tub), ribosomal protein L13 gene ( RPLI3) and 
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vacuolar-type ATPase gene ( V-ATPase) were cloned, and their expression levels in different 
developmental stages and different adult tissues were measured by qRT-PCR. Then the stabilities of these 
candidate genes were evaluated by a series of programs including ACt method, GeNorm, NormFinder, 
BestKeeper and an online program, RefFinder. Finally, the stabilities of selected reference genes were 
validated with an olfactory receptor co-receptor gene ( Orco) and pheromone binding protein 1 gene 
(PBP1). [Results] Analyzed by four programs including ACt method, GeNorm, NormFinder and 
BestKeeper, similar rankings of six candidate genes were obtained, among which RP49, RPL13 and 
GAPDH were the most stable genes in different developmental stages and different adult tissues, and ACT 
was ranked as the least stable gene despite of experimental conditions by all programs. Comprehensive 
analysis with RefFinder further showed that in different adult tissues, RPL13 was the most stable gene, 
followed by RP49, and in different developmental stages, RP49 was the most stable gene, followed by 
GAPDH. Additionally, the optimal number of reference genes was calculated by GeNorm as 2. Finally, 
the stabilities of selected reference genes were validated with Orco and PBP1 as target genes. The results 
indicated that when two pairs of genes, RPL13 and RP49, and RP49 and GAPDH, were respectively 
used as reference genes, Orco and PBP1 showed reliable expression patterns, which were consistent with 
life habits of C. punctiferalis and the results of the previous research. The expression patterns of Orco and 
PBPl, however, were irregular when ACT was used as the reference gene. [ Conclusion] In C. 
punctiferalis, two pairs of genes, RPL13 and RP49, and RP49 and GAPDH, are recommended as 
reference genes in different adult tissues and different developmental stages, respectively. 
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Dk AEE irt Conogethes punctiferails 广泛 分 布 于 东 
亚 南亚 、 东 南亚 和 澳大利亚 和 新 西 兰 地 区 ,其 幼虫 
为 害 多 种 果树 、 作 物 和 蔬菜 (Waterhouse, 1993) 。 在 
中 国 , 该 虫 的 寄主 植物 达 100 多 种 ,近年 来 为 害 玉米 
严重 ( 鹿 金 秋 等 , 2010) 。 桃 星 蜡 以 幼虫 钻 星 为 害 ， 
导致 对 该 虫 缺 乏 高 效 的 杀 虫 剂 ,同时 ,基于 性 信息 素 
的 性 诱 剂 产品 的 绿色 高 效 防 治 方式 ,已 经 在 桃 星 量 
的 防治 和 监测 上 得 到 广泛 应 用 ( 任 帅 和 郭 素 娟 ， 
2015), 

实时 荧光 定量 PCR ( quantitative real-time PCR, 
qRT-PCR) 以 其 可 靠 性 、 特 异性 和 高 灵敏 性 的 特征 ， 
在 过 去 的 几 十 年 中 被 广泛 用 于 基因 定量 分 析 ( 葛 星 
等 , 2013; 页 小 俭 等 , 2015; Jia et al., 2016) 。 在 
qRT-PCR 过 程 中 ,为 了 校准 样品 之 间 的 差异 ,确保 
定量 的 准确 性 ,需要 使 用 合适 的 内 参 基因 。 内 参 基 
因 往往 是 在 不 同 个 体 及 组 织 细胞 中 均 稳 定 表 达 的 看 
家 基因 ,这 些 基因 通常 参与 细胞 的 基本 生命 过 程 , 比 
如 微 管 蛋 白 (a-tub) 、 肌 动 蛋 白 (ACT) 甘油 醛 -3- 磷 
酸 脱氧 酶 (GAPDH) 和 18S rRNA 等 基因 (Daly et 
al., 1997; Thellin et al., 1999; Chen et al., 2005; Lu 
et al., 2013) 。 研 究 证 实 ,这 些 看 家 基因 的 转录 水 平 
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在 不 同 条 件 下 的 稳定 性 是 存在 差异 的 (Jonge， 
2007; Kozera and Rapacz, 2013) ,由 此 可 知 , 当 使 用 
内 参 基因 对 目标 基因 的 表达 进行 校准 时 ,内 参 基因 
自身 表达 的 稳定 性 必然 会 影响 定量 测定 的 准确 性 。 
因此 ,在 一 定 条 件 下 开展 基因 定量 研究 ,应 该 筛选 和 
验证 在 该 条 件 下 稳定 表达 的 内 参 基 因 。 另 外 ,在 内 
参 基因 的 个 数 方面 ,使 用 两 个 其 至 更 多 的 内 参 基因 ， 
能 够 获得 更 准确 的 定量 数据 ( Huang et al., 2015; 
Yang et al., 2015b; Ma et al., 2016) , 

近年 来 ,对 桃 星 旦 基因 方面 的 研究 逐渐 开展 ,其 
中 多 数 涉及 到 基因 的 定量 分 析 ( 葛 星 等 , 2013; 页 
小 俭 等 , 2015; Jia et al., 2016) ,然而 ,到 目前 为 止 
尚未 有 桃 蛙 虹 内 参 基 因 的 筛选 方面 的 报道 。 本 研究 
选取 桃 星 昌 不 同 发 育 时 期 和 成 虫 不 同 组织 两 种 条 
件 , 对 6 个 候选 基因 (4C7T，C4PD 厅 ，RP49 a-tub, 
RPL13 和 V-ATPase ) 的 转录 水 平 进 行 定 量 分 析 , 依 
据 内 参 基 因 筛 选 鉴定 的 常用 方法 ,使 用 ACt 法 ， 
GeNorm, NormFinder, BestKeeper 和 RefFinder 5 种 
软件 进行 稳定 性 分 析 ,最 后 ,对 筛选 出 的 候选 基因 的 
稳定 性 进行 验证 。 研 究 结 果 对 桃 星 蜡 的 基因 定量 研 
究 具 有 重要 意义 。 
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1 材料 与 方法 


1.1 试 虫 和 样品 

桃 星 蜡 室内 种 群 的 虫 源 采 自 重庆 市 合川 区 板栗 
园 ,在 室内 以 玉米 Zea mays 饲养 幼虫 (Xiao et al., 
2012). 。 饲 养 温 度 为 25 + 1%C ,相对 湿度 为 40% - 
60% , 光 周 期 为 13L:9D。 以 1596 的 蔗糖 水 供 成 虫 
补充 营养 。 收 集 桃 星 蜡 的 卵 (85 mg) 2 龄 (65 头 )、 
3 龄 (5 头 ) 5 龄 幼虫 (1 K) 分 别 混合 、 肾 (初期 ) 和 
成 虫 (3 HER) 虫 体 ( 艾 鹏 月 等 ,2014 ) 作为 不 同 发 育 
时 期 的 样品 。 收 集成 虫 头 (60 个 ) 、 胸 (15 个 )、 腹 (3 
A) KE CIS 对 ) 、 翅 (15 对 ) 和 触角 (45 对 ) 分 别 混 
合 。 将 所 有 样品 置 于 1.5 mL 无 RNA 酶 的 离心 管 
中 ,储存 于 -80% 条 件 下 用 于 RNA 提取 ,每 样品 准 
备 3 个 重复 。 
1.2 总 RNA 提取 和 cDNA 合成 

使 用 动物 组 织 总 RNA 提取 试剂 盒 (TIANCEN , 














中 国 ) 提 取 1.1 节 样 品 总 RNA。 进 行 DNA. 酶 处 理 
以 消除 总 RNA 中 的 基因 组 DNA 污染 ,最 后 溶 于 30 
pL ddH,0( AN RNase) 中 并 在 -80C 下 保存 。 提 
取 RNA 的 浓度 和 纯度 通过 BioTekEon 酶 标 仪 
(BioTek , 美国 ) 进行 测定 ,其 260/280 nm 的 RNA IR 
光度 比值 范围 为 1.8 ~2.1。 同 时 ,使 用 1% ZUR P 
凝 胶 电 泳 来 检测 RNA 完整 性 。 然 后 每 个 样品 取 
lpg RNA, 按照 PrimeScript™ RT t & 
( TakaRa&Clontech , 中 国 ) 要 求 的 步骤 合成 第 1 链 
cDNA 。 
1.3 ”内 参 基因 的 选择 和 克隆 

参考 其 他 昆虫 内 参 基因 (Lourenco et al., 2008; 
Lu et al., 2013; Zhai et al., 2014; Yang et al., 
2015a; Ma et al., 2016 ) ,选取 6 个 候选 基因 (AX 
1) 。 同 时 ,根据 桃 星 蜡 转录 组 的 序列 信息 (4C7 除 
外 ) ,并 通过 NCBI 比 对 其 他 蛾 类 中 同 源 基 因 的 相 
似 性 。 




















R1 候选 内 参 基因 信息 


Table 1 Candidate reference genes analyzed in this study 











GenBank 登录 号 




































































基因 基因 全 称 功能 
GenBank 
Genes Full gene name Function 
accession no. 
细胞 骨架 结构 蛋白 (调节 细胞 运动 ) 
ACT B- 肌 动 蛋白 B-Actin JX1190141 
Cytoskeletal structural protein (all cell; mediates cell motility) 
MESAN ERNIE SEE FEES PU 
Gin 甘油 醛 -3- 磷 酸 脱氧 酶 糖 酵 解 和 糖 异 生 作用 中 的 氧化 还 原 酶 KEERA 
Glyceraldehyde 3-phosphate dehydrogenase Oxidoreductase in glycolysis and gluconeogenesis 
RP49 核糖 体 蛋 白 49 Ribosomal protein 49 核糖 体 的 结构 成 分 Structural constituent of ribosome KX668533 
细胞 骨架 结构 蛋 E 
atub WEH Alpha tubulin RRS KX668534 
Cytoskeletal structural protein (flagella, median body, ventral disc ) 
崔 化 蛋白 质 合成 的 细胞 器 
RPL13 核糖 体 蛋 白 L13 Ribosomal protein L13 fe PE KX668535 
Organelles that catalyze protein synthesis 
H* 泵 入 细胞 管 泡 腔 内 或 泵 出 细胞 外 ,维持 细胞 内 正常 的 生理 
V-ATPase WWA! ATP 合 酶 Vacuolar-type ATPase 活动 KX668536 


Rotary molecular motors that generally function as proton pumps 


通过 Primer Premier 5 ( www. premierbiosoft. 
com) fll DNAMAN 8. 0 ( www. lynnon. com) 设计 克隆 
引物 ,同时 用 Primer 3 TF 385] 9] (4€ 2), PCR 扩 
增 为 25 pL (KA: 50 ng cDNA CUR EUR S 龄 期 幼 
Hh), 2.5 pL 10 x PCR 缓冲 液 ( 不 含 Mg * ), 62.5 
mmol/L Mg , 20 mmol/L dNTP, 0.5 U Taq ZZ jfi 

















(Takara, 中 国 ) ,引物 10 pmol/L, PCR 反应 条 件 : 
94°C 预 变性 3 min; 然后 94% 变性 30 s, 50% 退火 
30 s, 72% 延伸 1 min ,36 个 循环 ;最 后 72%C 5 min, 
最 后 保存 在 4% 。PCR 产物 经 1% BELT AM , 凝 
胶 提 取 试 剂 盒 进行 (OmegaBio-Tek，USA ) 纯化 , 然 
后 克隆 到 pGEM-T( Promega , R E] ) 中 进行 测序 。 
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R2 候选 内 参 基因 克隆 引物 信息 


Table 2 Information of cloning primers for candidate reference genes 





基因 引物 序列 (5 -37) 长 度 (bp) 
Genes Primer sequences Length 
ACT F: CACCCTGAAGTACCCCATCG R: TGAGGTAGTCGGTCAAGTC 776 
GAPDH F: ATCGGTCGCCTGGTGCT R: CGCCAATAAAGTCGGAAGA 822 
RP49 F: AAAGGTTTATCCGCCATCA R: TAGACGAGACGCCGTGAGC 263 
a-tub F: GGTCCGCACAGGCACAT R: GGGCGGCTGGTAGTTGAT 849 
RPL13 F: CCTGCCCGTGGACCTTT R: CCGCTTGACAGCCTTGA 305 
V-ATPase F: AGAATGAGGACAAGTTCGGTTAT R: GCCCACGATGGATACGG 1183 
1.4 qTR-PCR 检测 候选 内 参 基 因 行 综合 排名 。 
qRT-PCR 的 引物 信息 见 表 3。qRT-PCR 反应 体 ”1.6 内 参 基 因 稳 定性 验证 
系 为 20 uL: 2 x SYBR Green MasterMix ( Promega , 3 以 桃 性 旦 的 非典 型 气味 受 体 基因 Orco 和 性 信 


国 )10 pL, ddH,O 7 pL, 10 pmol/L 引物 ,cDNA 样 
m CB RC [n] Az A EST B e ake HAS [8] 2H 28.) 50. ngo 
每 个 样品 设置 2 个 技术 重复 和 3 个 生物 重复 。 扩 增 
步骤 :95Y 预 变性 2 min; 然后 进行 40 个 循环 :95% 
变性 15 s, 607C 退火 30 s; 最 后 在 607C 下 延伸 5 s。 
使 用 以 下 公式 计算 扩 增 效率 :E = [10^( - 1/slope ) 
-1] x 100% 。 获 得 具有 单 峰 的 炊 解 曲线 , 并 以 其 
确认 扩 增 后 每 个 引物 对 的 特异 性 以 及 是 否 具 有 引物 
ZRK, 
1.5 基因 表达 稳定 性 分 析 

通过 软件 GeNorm, NormFinder, 
RefFinder 以 及 相对 Act 法 ,根据 平均 循环 闽 值 (Ct 
值 ,3 个 生物 重复 ) 分 析 基 因 表 达 的 稳定 性 
( Vandesompele et al., 2002; Andersen et al., 2004; 
Pfaffl et al., 2004; Silver et al., 2006; Xie et al., 
2011), GeNorm 和 NormFinder 需要 输入 相对 定量 
值 ( 从 原始 Ct 值 转换 ) , 而 不 是 绝对 定量 循环 值 
(Cq) GeNorm 计算 内 参 基 因 表 达 稳 定性 的 M fr, 
根据 M 值 大 小 对 内 参 基因 的 稳定 性 进行 排名 ,同时 
GeNorm 还 可 以 计算 一 定 条 件 下 的 内 参 基 因 的 最 佳 
数目 (Vandesompele et al., 2002) 。NormFinder 则 通 
过 计算 不 同样 本 组 内 和 组 间 变 化 来 评估 基因 的 稳定 
性 。 另 外 ,根据 组 间 的 稳定 性 ,NormFinder 还 可 以 确 
定 两 个 基因 的 最 佳 组 合 。BestKeeper 根据 每 个 基因 
Cq 值 的 标准 偏差 和 变异 系数 对 各 个 基因 进行 排名 。 
而 ACt 方法 是 通过 两 两 比较 每 个 样品 中 基因 的 相对 
表达 量 以 鉴定 最 佳 基因 。 经 过 以 上 4 个 软件 分 析 之 
后 ,再 利用 在 线程 序 RefFinder (http; // fulxie. Ofees. 
us/? type = reference) ,对 所 有 候选 基因 的 稳定 性 进 






































BestKeeper , 


















































息 素 结合 蛋白 基因 PBPI 作为 目标 基因 ,分 别 以 软 
件 分 析出 的 最 稳定 /最 不 稳定 的 基因 作为 内 参 基因 ， 
使 用 2-~" 方 法 计算 Orco 和 PBP1 的 相对 表达 量 
(P«0.01, 邓肯 新 复 极 差 检 验 ) (Livak and 
Schmittgen, 2001) , 验证 候选 基因 的 稳定 性 。 











2 结果 


2.1 候选 基因 的 克隆 

通过 196 琼脂 糖 凝 胶 电 泳 检测 基因 的 扩 增 产 
物 ,经 纯化 和 克隆 后 ,通过 测序 获得 基因 序列 ,与 
NCBI 数据 库 比 对 ,确定 克隆 基因 的 正确 性 。 
2.2 引物 的 扩 增 效率 及 特异 性 评价 

qRT-PCR 结果 显示 ,每 个 候选 基因 的 溶解 曲线 
均 呈 现 单一 峰 , 表 明定 量 引 物 具有 良好 的 特异 性 。 
同时 ,基于 各 浓度 内 参 基因 的 cDNA 获得 的 标准 曲 
线 具 有 良好 的 决定 系数 (R >0.91) ,PCR 效率 也 在 
96. 8819% ~110.1300% 之 间 ( 表 3) ,表明 候选 基因 
的 定量 引物 设计 合理 ,能 够 用 于 相应 基因 的 定量 
测定 。 
2.3 候选 内 参 基因 的 表达 水 平分 析 

在 桃 蛙 坚 不 同 发 育 时 期 和 成 虫 不 同 组 织 中 两 个 
特定 条 件 下 ,每 个 候选 基因 的 Ct 值 如 图 1 所 示 。 所 
有 基因 的 Ct 值 均 高 于 10(14 ~22), 表 明 各 候选 基 
因 在 虫 体内 均 具 有 较 高 的 表达 水 平 。 同 时 各 个 基因 
的 表达 量 之 间 存 在 差异 ,其 中 4C7 A V-ATPase 的 Ct 
值 较 高 ,而 RP49 较 低 ,其 余 3 个 基因 居中 。 在 所 有 
的 6 个 基因 中 , RPLIS 表达 量变 化 最 小 , 其 次 是 
GAPDH 和 RP49 , {fii ACT 的 Ct 值 变化 最 大 。 
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表 3 候选 内 参 基因 qRT-PCR 引物 的 扩 增 效率 


Table 3 Amplification efficiency of primers used in qRT-PCR of candidate reference genes 

















基因 引物 序列 (5” -37) 扩 增 产物 长 度 (bp) 扩 增 效率 (% ) 线性 相关 系数 民 
Genes Primer sequences Amplification size Amplification efficiency Correlation coefficient 
F: CACCCTGAACTACCCCATCG 
ACT 374 96. 8819 0.984 
R: TGAGGTAGTCGGTCAAGTC 
F: ATGTCGTCGGAGTCAACCTT 
GAPDH 245 102. 5886 0.9912 
R; ACCAGTGGAAGCAGGAATGA 
F; GCGTAACTGGCGTAAACCC 
RP49 154 103.5765 0. 9836 
R: GCGGACGTTGTGTACTAGTA 
F: TCACTGGTAAGGAGGATGCG 
a-tub 238 100. 5972 0.9851 
R: TAGATCGCGAACTCCAGCTT 
F: GCTCCGAGTCTTCTGTCTGA 
RPL13 221 102. 7519 0.9902 
R: GGATGACTGCTGTGAAGGGA 
F: TGGTGGGAGAAATCATCCGT 
V-ATPase 186 103. 6685 0.993 
R: GCTCGTTGATGTCCTTCAGC 
F: ATGATGAACAAAGTGAAAGCC 
Orco 200 100. 0000 0.98335 
R: ACCCGGTACAAGTGACGGTACG 
F: ATGGAGCTGATGACGCCAC 
PBP!1 180 110. 1300 0.91784 
R: CCAAGATCTCTGCCAACACG 


基因 表达 量 (Ct 值 ) 
Gene expression level (Ct value) 





Co ACT 
ES a-iub 
EA GADPH 
EJ RP49 
RPL13 
E V-ATPase 


内 参 基因 Reference genes 














Fig. 1 


到 1 内 参 基 因 的 Ct 值 分 布 图 


Ct values distribution of reference genes 





上 图 表示 候选 内 参 基因 在 桃 性 蜡 不 同 发 育 时 期 及 成 虫 不 同 组 织 中 的 Ct 值 分 布 的 箱 型 图 ,样本 数 为 36。The figure shows the box diagram of Ct 


value distribution in different developmental stages and different adult tissues of Conogethes punctiferails ，and the number of samples was 36. 














2.4 不 同 条 件 下 内 参 基因 稳定 性 分 析 
2.4.1 ACt 和 GeNorm 分 析 : ACt 法 利用 基因 的 稳 
定性 与 ACt 值 的 标准 偏差 值 成 反比 的 原理 ,通过 对 
候选 基因 的 AC 值 的 标准 偏差 进行 两 两 比较 ,对 基 
性 进行 排序 。 结 果 显 示 , 在 成 虫 不 同 

织 中 , RPL13 的 稳定 性 最 高 , 其 次 是 RP49 和 
C4PDT。 同 时 ,最 不 稳定 的 基因 为 4C7( 图 2: A). 
而 在 不 同 发 育 时 期 中 ,RP49 则 为 最 稳定 的 基因 ,而 
ACT 同样 最 不 稳定 (图 2: B); 

与 ACt 法 相似 ,GeNorm 法 对 各 个 基因 稳定 性 生 
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RP49 ,在 不 同 发 育 时 期 则 是 RPA9 和 GAPDH, 而 最 
不 稳定 的 基因 则 仍然 是 4C7。 
2. 4. 2 NormFinder 分 析 和 BestKeeper 分 析 : 
NormFinder 通过 比较 样本 的 组 内 变异 和 组 间 变 异 来 
评估 基因 的 稳定 性 。NormFinder 分 析 表 明 , 在 不 同 
发 育 时 期 和 成 虫 不 同 组 织 两 个 条 件 下 ,最 稳定 均 为 
y RPL13 ,其 次 为 GAPDH, Tij ACT 则 为 最 不 稳定 基 
因 ( 图 3: A, B). 

Bestkeeper 则 根据 Ct 值 的 标准 差 来 分 析 基 因 的 
稳定 性 。 结 果 显 示 ,RP49 在 两 个 条 件 下 均 表 现 为 最 









































排序 也 是 依据 成 对 比较 ,所 以 GeNorm 的 排序 结果 
与 ACt 法 基本 一 致 。 如 图 2(C) 和 2(D) 所 示 , 在 成 
虫 不 同 组 织 中 最 稳定 的 基因 依次 是 RPL13 和 














稳定 的 基因 。 而 4C7 的 稳定 性 最 低 ,与 其 他 软件 的 
结果 相同 。 
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Fig. 2 Stability analysis of candidate reference genes in different adult tissues (A and C) and different 























developmental stages ( B and D) of Conogethes punctiferails 


A 和 B 采 用 ACt 法 ,C H D RH GeNorm 软件 分 析 。ACt method was used for A and B, and GeNorm program for C and D. 
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Fig. 3 Stability analysis of candidate reference genes in different adult tissues (A and C) and different 
























































developmental stages ( B and D) of Conogethes punctiferails 





A All B RJ NormFinder 软件 分 析 ,C 和 D 采用 BestKeeper 软件 分 析 。NormFinder program was used for A and B, and BestKeeper program for C 


and D. 


2.4.3 RefFinder 分 析 : 经 过 以 上 4 个 软件 的 独立 。 其 次 为 RP49 和 GAPDH ,最 不 稳定 的 为 4CT( 图 4: 


分 析 之 后 ,再 利用 在 线程 序 RefFinder 综合 评估 内 参 A) 。 对 于 不 同 发 育 时 期 ,最 稳定 的 基因 为 RP49， 
基因 的 稳定 性 并 进行 排名 (Xie et al., 2011) 。 结 果 ”次 为 G4PDH 和 RPL13 ,ACT 则 依然 为 最 不 稳定 的 
显示 ,在 成 虫 不 同 组 织 中 , 最 稳定 的 基因 为 RPL13， 因 ( 图 4: B) 
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图 4 RefFinder 软件 分 析 桃 性 量 成 虫 不 同 组 织 (A) 和 不 同 发 育 时 期 (B) 候 选 内 参 基 因 的 稳定 性 
Fig. 4 Stability analysis of candidate reference genes in different adult tissues (A) and different developmental 
stages (B) of Conogethes punctiferails by using RefFinder analysis 
2.5 ”内参 基因 最 佳 数目 的 确定 于 0.15 ,说 明 内 参 基因 的 最 佳 数 目 为 2, 在 以 上 两 种 








利用 GeNorm 软件 通过 判断 配对 变异 值 ( 低 于 
0. 15) 来 确定 内 参 基因 最 佳 数目 , 即 当 Vn/(n +1) 
«0.15 时 ,说 明 该 条 件 下 最 佳 内 参 基因 的 数目 为 n 
个 。 软 件 分 析 结 果 如 图 5 所 示 , 在 不 同 发 育 时 期 和 
成 虫 不 同 组 织 条 件 下 ,只 有 V2/3 的 配对 变异 值 低 

















条 件 下 进行 基因 定量 工作 时 ,最 好 使 用 2 个 内 参 基 
因 。 因 此 ,结合 以 上 内 参 基因 稳定 性 分 析 的 结果 可 
知 ,在 成 虫 不 同 组 织 条 件 下 ,可 以 使 用 RP49 和 
RPL13 作为 内 参 基 因 ,而 在 不 同 发 育 时 期 中 则 可 以 
使 用 RP49 和 GAPDH 作为 内 参 基因 。 
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Fig. 5 Evaluation of the optimal number of reference genes for normalization in different adult tissues 


and developmental stages of Conogethes punctiferalis 


虚线 表示 配对 变异 值 为 0. 15。The dashed line indicates that the pairwise variation is 0. 15. 


2.6 内参 基因 稳定 性 的 验证 

在 对 候选 内 参 基因 的 稳定 性 进行 软件 分 析 之 
后 ,需要 对 分 析 结 果 进 行 验证 。 在 本 实验 中 ,首先 选 
择 在 昆虫 中 广泛 分 布 的 嗅觉 关键 基因 Orco 作为 目 
的 基因 ,以 候选 基因 (不 同 发 育 时 期 : GAPDH 和 
RP49 ;成 虫 组 织 : RP49 和 RPLA3 ,以 及 最 不 稳定 的 
基因 4C7) 作为 内 参 基 因 ,测定 目的 基因 在 不 同 发 育 
时 期 和 不 同 组 织 两 种 条 件 下 的 表达 量 。 结 果 如 图 6 



































(A) 所 示 , 以 RP49 和 RPL13 为 内 参 基 因 时 ,目的 基 
Orco 表达 量 在 成 虫 的 6 个 不 同 组 织 之 间 均 无 显 
著 性 差异 ,在 触角 中 的 表达 量 远 高 于 其 他 组 织 ; 以 
ACT 为 内 参 基因 时 , Orco 表达 量 在 头 部 及 翅 部 均 分 
别 显 著 高 于 以 RP49 和 RPL13 为 内 参 基 因 时 的 表达 
量 。 相 似 的 结果 也 出 现在 不 同 发 育 时 期 的 条 件 下 
(图 6: B), DJ GAPDH, RP49 JJ ACT 为 内 参 基因 
时 ，Orco 基因 在 成 虫 期 的 表达 量 均 为 其 他 发 育 时 期 
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R6 Orco 在 桃 蛙 坚 成 虫 不 同 组 织 中 (A) 和 不 同 发 育 时 期 (B) 的 表达 水 平 





Fig. 6 Expression level of Orco in different adult tissues (A) and different developmental stages (B) of Conogethes punctiferalis 


A 和 B 图 分 别 以 RP49, RPLI3 和 ACT 及 GAPDH, RP49 Fil ACT HWER 
































因 。 图 中 数据 为 平均 值 + 标准 差 。 柱 上 星 号 示 同 一 组 织 或 发 育 时 期 


不 同 内 参 基 因 间 差异 显著 (已 <0.01, 邓肯 新 复 极 差 检验 ) 。RP49, RPLI3 and ACT (A) and GAPDH, RP49 and ACT (B) were used as reference 


genes in Figs. A and B, respectively. Data in the figure are mean + SD. Asterisk above bars means significant difference ( P <0.01, Duncan' s test) 





between different internal reference genes in the same tissue or the same developmental stage. 


(Bp ,2 龄 幼虫 ,3 龄 幼虫 ,3 龄 幼虫 和 晴 ) 的 2 倍 以 
上 ;以 GAPDH 和 RP49 为 内 参 基 因 时 , Orco 的 表达 
量 在 每 个 发 育 阶段 之 间 均 无 显著 差异 ,以 ACT 为 内 
BIN Orco 表达 量 而 在 卵 及 肾 期 均 显 著 高 于 以 
GAPDH 和 RP49 为 内 参 基因 时 的 表达 量 。 

为 客观 反映 最 适 内 参 基 因 的 适用 性 ,以 PBPI 
为 目的 基因 ,进一步 检测 了 RPA9, RPLI3, GAPDH 
All ACT 4 个 内 参 基因 在 桃 蛙 螟 成 虫 不 同 组织 中 和 
不 同 发 育 时 期 的 表达 量 。 结 果 如 图 7(A) 所 示 , 以 





RP49 和 RPLI3 为 内 参 基因 测 得 的 PBP1 表达 量 在 
成 虫 的 6 个 不 同 组 织 之 间 均 无 显著 差异 ; 以 ACT 为 
内 参 基 因 时 ,PBP1 在 头 部 和 腹部 的 表达 量 分 别 显 
著 高 于 以 RP49 和 RPL13 为 内 参 基 因 时 的 表达 量 。 
由 图 7(B) 可 见 ,在 以 GAPDH 和 RP49 为 内 参 基因 
测 得 的 PBP1 表达 量 在 卵 .2 龄 幼虫 3 龄 幼虫 .5 龄 
幼虫 和 肾 之 间 均 无 显著 差异 ,以 ACT 为 内 参 基因 
时 ,PBP1 在 卵 及 3 龄 幼虫 期 的 表达 量 均 显 著 高 于 
以 GAPDH 和 RP49 为 内 参 基 因 的 表达 量 。 
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Fig. 7 Expression level of PBP1 in different adult tissues (A) and different developmental stages ( B) of Conogethes punctiferalis 


A 和 B 图 分 别 以 RP49, RPL13 Al ACT 及 GAPDH, RP49 Fil ACT HWER 

















因 。 柱 上 星 号 示 同 一 组 织 或 发 育 时 期 不 同 内 参 基因 间 差 异 显著 (P < 





0.01, 邓肯 新 复 极 差 检 验 )。RP49,， RPLI3 and ACT (A) and GAPDH, RP49 and ACT (B) were used as reference genes in Figs. A and B, 


respectively. Data in the figure are mean + SD. Asterisk above bars means significant difference ( P «0. 01, Duncan’s test) between different internal 


reference genes in the same tissue or the same developmental stage. 


3 讨论 





qRT-PCR 是 基因 表达 分 析 研 究 中 的 重要 技术 ， 
该 技术 采用 相对 定量 的 方法 对 基因 的 表达 量 进行 检 
测 (Guan and Yang, 2008), 即 以 内 参 基因 作为 标准 
校准 目的 基因 的 表达 量 , 具 有 灵敏 度 高 ,重复 性 好 的 














特点 , 目前 已 成 为 基因 表达 研究 中 检测 或 比较 
mRNA 水 平 的 标准 技术 。 目 前 在 多 种 昆虫 中 都 进行 
itp BS FEB AY Si 2 , HH AE AE Bombyx mori ( Peng et 
al., 2012) , 18 KEL Nilaparvata lugens ( Yuan et al., 
2014) , ZR JUL X. Tribolium castaneum ( Xl] & 18 |, 
2014) 、 小菜 峨 Plutella Xylostella( 符 伟 等 , 2012) 二 
BEN i Tetranychus urticae ( Morales et al., 2016) FIJA 
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SCE Bactrocera Cucurbitae ( 周 世 豪 等 , 2016 ) 。 本 研 测 得 的 Orco 及 PBPl 表达 量 , 具 有 一 致 的 表达 趋 


究 基 于 桃 星 蜡 转 录 组 的 测序 数据 (Xiao et al., 
2016) ,并 参考 其 他 昆虫 常用 的 候选 内 参 基因 ,通过 
软件 分 析 基 因 表 达 稳 定性 ,在 不 同 发 育 时 期 和 成 虫 
不 同 身 体 部 位 两 种 条 件 下 , 筛 选 出 桃 星 曝 的 内 参 基 
o 

候选 内 参 基 因 的 稳定 性 分 析 , 常 采用 系列 软件 
依据 相关 参数 分 析 各 个 基因 的 稳定 性 ,并 根据 稳定 
性 高 低 进 行 排序 。 本 研究 首先 使 用 4 个 软件 进行 独 
立 排序 ,发 现 RP49, GAPDH 和 RPL13 是 6 个 候选 
基因 中 最 稳定 的 3 个 基因 , 只 是 由 于 不 同 软件 的 计 
算 原 理 不 相同 ,给 出 的 3 个 基因 的 稳定 性 高 低 略 有 
差异 (图 2 和 3)。 进 一 步 通过 在 线程 序 RefFinder 
给 出 了 候选 基因 的 综合 排名 :在 不 同 发 育 时 期 ,最 稳 
定 的 基因 为 RP49 ,其 次 为 C4PDH。 而 在 成 虫 不 同 
身体 部 位 条 件 下 ,最 稳定 的 基因 是 RPL, 其 次 为 
RP49( 图 4)。 在 其 他 昆虫 或 螨 中 ,RPL13,，RP49 和 
GAPDH 的 表达 量 通 常 比较 稳定 ( Lourenço et al., 
2008; Van Hiel et al., 2009; Yang et al., 2015a) , 
另外 ,对 于 稳定 性 最 差 的 基因 ,4 个 软件 的 单独 排序 
结果 与 在 线程 序 RefFinder 的 综合 排序 结果 完全 相 
le] , 均 为 4C7。 该 基因 在 樱桃 果 晶 Drosophila suzukii 
TI Aq HT ios Tetranychus cinnabarinus mite 的 内 参 基 
因 筛 选中 ,也 被 评估 为 最 不 稳定 的 基因 (Sun et al., 
2010; Zhai et al., 2014) 。 当 然 ,在 有 的 物种 中 ACT 
也 被 评估 为 最 佳 内 参 基 因 (Lu et al., 2013; Shang et 
al., 2015) 。 这 些 结果 进一步 表明 ,内参 基 因 的 表达 
稳定 性 因 不 同 物种 而 不 同 ,同时 也 受到 实验 条 件 的 
影响 。 

在 分 析 了 候选 基因 表达 的 稳定 性 之 后 ,进一步 
使 用 GeNorm 软件 计算 出 桃 蛙 晨 内 参 基因 的 最 佳 个 
数 为 2 个 。 在 基因 表达 研究 中 ,通常 使 用 单一 的 内 
参 基因 ( Deng et al., 2016; Liu et al., 2016)。 然 而 ， 
近年 来 ,多 个 内 参 基 因 被 广泛 应 用 ,由 此 能 提高 对 目 
的 基因 定量 的 准确 性 (Vandesompele et al., 2002) 。 
因此 ,在 桃 星 蜡 的 基因 定量 研究 中 ,为 提高 定量 的 准 
确 性 ,在 不 同 发 育 时 期 和 成 虫 不 同 组 织 的 条 件 下 ,也 
推荐 使 用 2 个 内 参 基 因 。 因 此 ,建议 将 RPL13 和 
RP49 作为 成 虫 不 同 身体 组 织 条 件 下 的 内 参 基因 ,而 
将 RP49 和 CAPDH 作为 不 同 发 育 时 期 条 件 下 的 内 
参 基因 。 

最 后 ,以 重要 的 嗅觉 受 体 基因 Orco 和 PBP1 验 
证 了 两 对 基因 RPL13 和 RP49 以 及 RP49 和 GAPDH 
的 稳定 性 。 结 果 发 现 ,分 别 以 两 对 内 参 基 因为 基准 

























































































势 , 即 在 触角 和 成 虫 阶段 高 表达 ,在 其 余 组 织 中 和 发 
育 时 期 的 表达 量 相对 较 低 (图 6: A, B) 。 在 蛾 类 昆 
虫 的 嗅觉 中 ,成虫 阶段 以 触角 为 主 的 嗅觉 行为 具有 
非常 重要 的 作用 (Stengl and Funk, 2013; Carraher et 
al., 2015) ,在 本 实验 中 ,采用 RPL13 和 RPA9 以 及 
RP49 和 GAPDH 两 对 内 参 基因 对 Orco 和 PBP1 的 定 
量 结果 符合 蛾 类 昆虫 的 嗅觉 规律 。 而 男 一 方面 , 当 
使 用 4C7 作为 内 参 时 ,Orco 和 PBPL 的 定量 结果 在 
两 个 条 件 下 的 均 不 稳定 ,说 明 4C7 确实 不 适合 作为 
桃 星 蜡 的 内 参 基 因 。 

ARE FE WA AE HE dE 5 龄 幼虫 中 克隆 了 ACT, 
GAPDH , RP49 , a-tub, RPL13 和 V-ATPase 6 4- 
基因 ,通过 ACt YE, GeNom, 
NormFinder 和 RefFinder 软件 ,分析 了 各 个 基因 的 稳 
定性 。 结 果 发 现在 不 同 发 育 时 期 ,最 稳定 的 基因 为 
RP49 ,其 次 为 G4PDH。 而 在 成 虫 不 同 组 织 中 ,最 稳 
定 的 基因 是 RPL13 ,其 次 为 RP49。 在 两 种 条 件 下 进 
行 桃 星 蜡 的 基因 定量 研究 时 ,为 得 到 更 准确 的 定量 
结果 ,建议 使 用 2 个 内 参 基 因 。 在 不 同 发 育 时 期 条 
件 下 建议 使 用 RP49 和 C4PD 末 ,在 成 虫 不 同 身体 部 
位 条 件 下 建议 使 用 RPLI3 和 RP49。 本 研究 结果 有 
利于 进一步 揭示 桃 蛙 蜡 嗅觉 相关 的 分 子 机 理 , 即 定 
量 分 析 嗅 觉 基 因 的 mRNA 表达 水 平 ,从 而 为 后 续 
RNA 干扰 等 相关 研究 奠定 基础 。 
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